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Abstract : An analysis is performed of the liquid-gas phase transition of nuclear matter 
obtained from different versions of scalar derivate coupling suggested by Zimanyi and 
Moszkowski (ZM) and the results are compared with those obtained from the Walecka 
model. We present the phase diagram for the models and one of them, the ZM3 model, 
has the lowest critical temperature = 13.6 MeV with the lowest critical density pc = 
0.037 fm~^ and pressure Pc = 0.157 MeV fm-3. These results are in accord with recent 
observations from energetic heavy-ion collisions, which suggest a small liquid-gas phase 
region. 
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1 Introduction 



Nowadays, the study of the hquid-gas phase transition, which may occur in the warm 
and dilute matter produced in energetic heavy-ion colhsions, is one of the interesting prob- 
lems in nuclear physics [|l|. This idea of that nuclear systems may show up a critical be- 
havior has initiated more than ten years ago with the observation by the Purdue-Fermilab 
group of asymptotic fragment charge distributions exhibiting a power law 0. This in- 
terest increased recently with the attempt by the EOS Collaboration to extract critical 
exponents of fragmenting nuclear systems produced in the collision of 1 GeV/nucleon Au 
nuclei with a carbon target [Q], and with the extraction by ALADIN/LAND Collaboration 
of a caloric curve resulting from the fragmentation of the quasiprojectile formed in the 
collision Au + Au at 600 MeV/nucleon exhibiting a behavior expected for a first order 
liquid-gas phase transition 

At the time where the search for signals of liquid-gas nuclear phase transition are 
taking place, it is important to have ready the theoretical phase-transition predictions 
for a broad class of different hadronic models. The main ingredient in this analysis 
is the nuclear matter equation of state (EOS) at finite temperature. The success of 
relativistic mean-field theories describing cold nuclear matter and bulk nuclear properties 
throughout the periodic table, suggests the use of a relativistic mean-field EOS. Moreover, 
the mean-field approximation is known to be thermodinamicaly consistent: the relevant 
thermodynamic identities are all satisfied 0, p. 

Recently variants of Zimanyi and Moszkowski (ZM) model were implemented and 
applied by us to dense and cold nuclear matter [§], [|]. The usual ZM model, also referred 
to in the literature as the Derivative Scalar Coupling model (DSC), consists of derivative 
coupling between nucleons and scalar mesons a. The model has been extended to include 
a non-linear interaction between the nucleon and the vector meson uo. Two types of 
this interaction were employed and the resulting models were denoted ZM2 and ZM3. 
These models were designed to cure the defects of the Walecka model [0, namely the 
low effective nucleon mass and the large incompressibility of nuclear matter. Each one of 
them is very simple since they have only two free parameters, the scalar (vector) coupling 
constants (C^), adjusted to reproduce the binding energy {Eb) of the nuclear matter 
at p = po- The degrees of freedom are baryon fields [ip), scalar meson fields (a), and 
vector meson fields {uj). 

In all ZM models, there are non-linear interaction terms which in a approximate 
way incorporate the effect of many-body forces. After an appropriate rescaling of the La- 
grangians, these models can be understood as generalizations of the Walecka Model where 
the scalar and vector meson couplings become effectively density-dependent |TT| . This fact 



underlies the recent approach, known as relativistic density-dependent Hartree-Fock |12 
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|T3[, which describes finite nuclei and nuclear matter saturation properties using cou- 
pling constants that are fitted, at each density value, to the relativistic Brueckner-Hartree- 
Fock self-energy terms. The good agreement obtained for the ground state properties of 
spherical nuclei lends support to such density-dependent coupling constants. Recently, it 
was shown that chiral symmetry restoration requires the meson-nucleon coupling to be 
density dependent ||15|. 

The aim of this paper is to extend our study to include temperature effects, and to 
perform an analysis of the liquid-gas phase transition of the warm nuclear matter obtained 
on these three ZM models and compare to the linear Walecka (W) model [|l6l. We present 
the effective nucleon mass, energy per nucleon, pressure, and entropy density as a function 
of the baryonic density at different temperatures. We show the isotherms, construct the 
phase diagram with the phase coexistence boundary, and present the critical and flash 
temperatures for the models. The usual ZM model has already been applied to investigate 
some thermodynamic properties of nuclear matter []T7| and recently, the modifled versions 
have been used in order to study the density and temperature dependence of hadron 
masses 
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The outline of the paper is as follows: in the next section we present the EOS at 
finite temperature. Section 3 includes our results and discussion of the thermodynamic 
properties of nuclear matter. Finally, we summarize. 



2 The nuclear matter EOS at finite temperature 

Since the models we are dealing with were discussed in detail in references [0, 0, 0, 
here we will only present the Lagrangian obtained after rescaling the nucleon field as 
ijj m*^^'^i/j for all ZM models and making the rescaling tu^ Tn*uj^j, for ZM2 and ZM3 
models: 

- ^lj{M-m*'g^a)tP + ^id,crd''a-mla'), (1) 

where a and /3 have the following values for the different models, W : a = 0,(3 = 0; 
ZM : a = 0, /5 = 1; ZM2 : a = 1, /5 = 1; ZM3 : a = 2, /5 = 1; and m* = (1 + g^a/M)'^ 
in all three cases, M is the bare nucleon mass and F^jy = d^Uu — dyOJ^j,. 

When the meson fields are replaced by the constant classical fields cxo and uJq we arrive 
at the mean-field approximation, with the equation of motion for the nucleon: 

[^7^5^ - (M - m*'g,a) - m*" g^-f,u^]^ = 0, (2) 
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where the effective nucleoli mass M* is given by M* = M — m*^gcr(J. In the case of ZM 
models where /3 = 1 we can identify m* = M* /M = (1 + gf^a / M)~^ . 

The expression for the energy density and pressure at a given temperature T can be 
found as usual by the average of the energy-momentum tensor, 

P = 3,m*"p' - ^ C-^y + It^. I + (4) 



2M2 2Cl \ m*^ J 3 {2Txf J E*{k) 

Thus we obtain the entropy density: 



1 
T 



" ^* V + 7^ jd^kE*{k) {Uk + fik) 



M2 ^ (27r) 



+ — ;r- / CL k ^ ,, , (Uk + rik) —., (5) 

where 7 is the degeneracy factor (7 = 4 for nuclear matter and 7 = 2 for pure neutron 
matter ), and stand for the Fermi-Dirac distribution for baryons and antibaryons 
respectively, with arguments {E* — v)/T , E*{k) is given by E*{k) = (A;2 + M*2)2 . An 
effective chemical potential which preserves the number of baryons and antibaryons in 
the ensemble is defined hj u = ^ — V, with fi is the thermodynamical chemical potential. 
We have introduced = glM'^/ml and = g^M'^/m'^. 

The effective mass is obtained explicitly through the minimization of S with respect 
to m* and must satisfied the self-consistent equation, 

7C2 ^3^+1 f x'^dx aClCl ^o.+20 2 



where we have used the dimensionless variable x = -g. 

The energy density can be fitted to the nuclear matter ground state energy and sat- 
uration density po at zero temperature to obtain the different coupling constants for the 
models. They are presented in table 1 together with the nuclear matter incompressibility 
that at T = is given by: 

K = ^Pog^(-^\p=Po = ^Pog^\p=Po- (7) 

To compute the thermodynamic functions, one first chooses T and u. The self- 
consistency condition in eq.(^) is then solved to determine M* ( note that there may 
be several solutions for fixed T and u ). These solutions specify the distribution functions 
Uk and Uk, and the remaining integrals in Eqs. (^), (0) and (|^) can then be evaluated 
directly. 
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3 Results and Discussion 



In Fig. 1 we show M* as a function of T at zero density. In this regime, the vector 
field proportional to p vanishes, and so the three ZM models differ only in having different 
values of their scalar coupling constants C^. The ZM and the Walecka model coincide 
in the lower temperature region T < 120 MeV and ZM3 model stay together up to 
T ~ 160 MeV. However, at higher temperature the models separate quite clearly, with 
the effective nucleon mass in the ZM models dropping more slowly than that in the 
Walecka Model. This means that the sigma field (the source for the scalar density), 
increases more slowly with temperature in ZM models because of the inclusion of non 
linear interactions which are absent in the Walecka Model. As a result, the attraction 
is stronger in the Walecka model favouring the formation of nucleon-antinucleon pairs at 
high temperature. Moreover, none of the proposed ZM models is able to present a first 
order phase transition at p= , T ^0. This is in contrast to the Walecka model, which 
has such a phase transition at T ~ 185 MeV [p!9| . 

In Fig. 2 we show the behavior of the effective nucleon mass with density at different 
temperatures for all the models. For low temperatures the results are not so different from 
those obtained at zero temperature, showing that in this regime the density dependence 
is more important than the temperature dependence. As temperature is raised, M* first 
increases and then decreases more slowly for ZM models than Walecka model at T = 200 
MeV. Within ZM models, this decrease is more pronounced in ZM3, but is even smaller 
compared to Walecka model where the effective mass goes down very fast. In short, 
the effect of the temperature on the effective nucleon mass in the ZM models is not so 
pronounced as in the case of Walecka Model, and can be seen only for densities below the 
normal density. 

We present the energy per nucleon as a function of the density at various temperatures 
in Fig. 3 . As the temperature increases the nuclear matter becomes less bound and the 
the saturation curve around the equilibrium point in the ZM models is flatter than that 
in the Walecka model. This indicates that the nuclear matter EOS in ZM models is softer 
compared to the obtained in the Walecka Model, even at finite temperature. We can also 
conclude that the incompressibility of nuclear matter decreases when the temperature 
increases. This can be seen more clearly in Fig. 4 where we show the pressure-density 
isotherms of nuclear matter at different temperatures. Since the incompressibility K is 
related to dp/ dp (calculated at the equilibrium point where the pressure vanishes), we see 
directly that when the temperature increases K decreases, and among the ZM models, 
the ZM3 model always gives the softest EOS for a fixed temperature. 

The isotherms exhibit a typical Van der Waals like interaction where a liquid and 
gaseous phases coexist, with an unphysical region in the mi die of each isotherm that gets 
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smaller as the temperature increases. For very small temperatures the isotherms manifest 
the following behavior: for very low density the pressure increases with temperature as for 
as a ideal gas, p ~ pkbT. It decreases subsequently because of the attractive interaction 
of the sigma field, and finally increases as a consequence of the repulsion coming from the 
vector meson which dominates at high density. When temperature increases, the term 
pkhT becomes more important and the local minimum in the pressure is less pronounced 
and disappears when the temperature is equal to the critical Tc. At this temperature, the 
unphysical region disappears and an inflection point appears in the isotherm, as we show 
in the Fig.4 for each model. The p-p isotherms in the ZM models have a shallower and 
more even valley than the corresponding ones in the Walecka Model, and this is more 
noticeable in ZM3 model. In table 2 we list the critical temperature Tc, density pc and 
pressure pc given by the ZM and Walecka models. The ZM3 model presents the lowest 
Tc = 13.6 MeV, density pc = 0.037 fm~^ and pressure Pc = 0.157 MeV fm~^. 

The phase coexistence boundary is determined by Gibbs's criteria, namely, that the 
liquid and gas phases have equal temperatures (thermal equilibrium), chemical potentials 
(chemical equilibrium), and pressures (hydrostatic equilibrium). We present in Fig.5 the 
phase diagram Txp of the models. Below the coexistence curve of each one, the equilibrium 
state is a mixture of gas and liquid. This region is bigger in the Walecka model. In fact, 
if we include nonlinear terms in this model, this region becames smaller and the critical 
temperature goes down to Tc = 14.2 MeV 0. The ZM3 model, where the non-linearity 
of the coupling between the vector field to the nucleon is strongest, presents the smallest 
phase coexistence region comparing to the other models. 

As we have already pointed out, the nuclear matter incompressibility K decreases when 
the temperature increases. So, we will have a temperature where the incompressibility 
K calculated at the equilibrium point vanishes. This temperature is known as the flash 
temperature T = T^, = p(p/,T;) = 0. It represents the highest temperature 

at which a self-bound system can exist in hydrostatic equilibrium (p = 0). Above this 
temperature the warm nuclear matter is unbound and starts expanding. We present in 
Fig. 6 the pressure as a function of baryon density at the flash temperature for the 
models. This temperature in MeV is 14.1, 12.9, 12.2 and 11.0 for the Walecka, ZM, ZM2 
and ZM3 models respectively. Again, the ZM3 model has the smallest flash temperature. 
As expected, all of these temperatures are lower than the critical ones, because, as Fig.4 
shows, at the critical temperature the pressure is already positive and the system is 
expanding. 

Finally, we present in Fig. 7 the entropy density as a function of the density at different 
temperatures. For high temperatures (T = 200 MeV), we see an increase in the entropy 
density with the density for all the models. This happens even at very low densities, 
and manifests what we have already pointed out when we have discussed the behavior of 
the effective nucleon mass with the temperature at zero density. This decrease of M* or 
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increase of the entropy density with increasing temperature, which is more pronounced 
in the Walecka and ZM3 models, resembles a phase transition. At high temperature and 
low density the system becomes a dilute gas of baryons in a sea of baryons-antibaryons. 

In summary, we have presented the thermodynamic properties of nuclear matter in 
three different versions of the ZM model. We have shown how the effective nucleon mass 
M*, energy per nucleon, pressure, and entropy behaves as a function of the density for 
different temperatures. As in the zero temperature case, all the ZM models give a softer 
EOS of nuclear matter at finite temperature than the Walecka model. Among the three 
ZM models, ZM3 is the softest. Unlike the Walecka Model the ZM models do not exhibit 
a phase transition for finite temperature at zero density. We studied the liquid-gas phase 
transition, and concluded that the ZM3 model presents the smallest phase coexistence 
region with the lowest critical temperature, density and pressure. The incompressiblity 
decreases with the increasing temperature, and vanishes when T = Tfiash- Again, the 
ZM3 model has the smallest flash temperature. The experimental results investigating 
this warm and dilute matter produced in energetic heavy-ion collisions suggest a small 
hquid-gas phase region and a low critical temperature 0] , . This suggests that the good 
description of nuclear matter properties obtained in the ZM3 model at zero temperature 
remains, even at finite temperature, and makes the ZM3 model the most suitable of all 
models used. 

One of us (MM) would like to thank Dr. D. R. Phillips for careful reading of the 
manuscript, the Maryland TQHN Group supported by the U.S. Depto. of Energy for 
their hospitality during his extended stay, and gratefully the CAPES of the Brazilian 
Governement for the financial support which made this visit possible. 
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4 Figures and Captions 

Fig. 1: Baryon effective mass in nuclear matter as a function of the temperature at p = 0. 



Fig. 2: Baryon effective mass M* as a function of tlie baryon density at different temper- 
atures for the Walecka model (W) and Zimanyi-Moszkowski models (ZM, ZM2, ZM3). 

Fig. 3: Proper energy/baryon as a function of baryon density at different temperatures 
for the Walecka model (W) and Zimanyi-Moszkowski models (ZM, ZM2, ZM3). 



Fig. 4: Pressure as a function of baryon density at different temperatures for the Walecka 
model (W) and Zimanyi-Moszkowski models (ZM, ZM2, ZM3). 



Fig. 5: Temperature as a function of the baryon density (phase diagram) for the Walecka 
model (W) and Zimanyi-Moszkowski models (ZM, ZM2, ZM3). 



Fig. 6: Pressure as a function of baryon density at flash temperature (T/) for the Walecka 
model (W) and Zimanyi-Moszkowski models (ZM, ZM2, ZM3). 



Fig. 7: Entropy density as a function of baryon density at different temperatures for 
the Walecka model (W) and Zimanyi-Moszkowski models (ZM, ZM2, ZM3). 
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Table 1: Coupling constants and C^; binding energy Eb (MeV) at equilibrium density 
Po{fm~^), m* and the incompressibility K for the indicated models. 



models 




^2 


Eb 


Po 


m* 


K 


W 


357.4 


273.8 


-15.75 


0.148 


0.54 


550.82 


ZM 


169.2 


59.1 


-15.90 


0.160 


0.85 


224.71 


ZM2 


219.3 


100.5 


-15.77 


0.152 


0.82 


198.32 


ZM3 


443.3 


305.5 


-15.76 


0.149 


0.72 


155.74 



Table 2: Values for the critical temperature Tc and the effective mass M* in MeV, critical 
density pc in fm~^ and pressure Pc in MeV/ fw? for the indicated models. 



Models 


T 

c 


Pc 


Pc 


m; 


Walecka 


18.3 


0.0650 


0.4300 


760 


ZM 


16.5 


0.0698 


0.2570 


861 


ZM2 


15.5 


0.0364 


0.2106 


881 


ZM3 


13.6 


0.0354 


0.1571 


831 
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